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Abstract 

(N 

ff") ' We study the decay process of the Z boson to a photon and a graviton. The most general 

form of the on-shell amplitude, subject to the constraints due to the conservation of the 
electromagnetic current and the energy-momentum tensor, is determined. The amplitude 
is expressed in terms of three form factors, two of which are CP-odd while one is CP-even. 
The latter, which is the only non-zero form factor at the one-loop level, is computed in the 
standard model and the decay rate is determined. 
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l> 1 Introduction 

As is well known pQ, the Z-boson cannot decay into two photons. Therefore, the simplest decay 
^ \ of the Z-boson to two massless bosons is through the channel 

Z{p) _ 7 (fc) + Q{q) , (1.1) 

where 7 denotes the photon and Q the graviton. In this work we consider this process and 
calculate the decay rate, assuming the standard electroweak interactions and the canonical 
gravitational coupling of the standard particles. 

We work with the linear theory of gravity, which means that we write the space-time metric 
in the form 



9\ P = r]\ P + 2Kh Xp (1.2) 

where h\ p is identified with the graviton field, and the gravitational couplings in the Lagrangian 
are expanded up to the linear order in k. Furthermore, the constant k is defined in terms of 
Newton's gravitational constant by 

K = y/8nG, (1.3) 

which is such that the field h\ p has the properly normalized kinetic energy term in the La- 
grangian. This point of view for treating processes involving the gravitational and Standard 
Model interactions is the same as that employed in some recent works for the calculation of 
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quantum gravity amplitudes [U |3] , in which General Relativity is treated as an effective field 
theory for energies below the Planck scale [31 El IH] - 

The process is a clean one as far the kinematics is concerned: it involves the photon, with 
an energy equal to half the Z-boson mass. However, the rate is not sizable, as a simple order- 
of- magnitude estimate of the decay rate readily shows. The gravitational couplings will involve 
the factor k in the amplitude, and therefore G in the rate. The Z and the photon couplings 
provide a factor a 2 in the rate. Dimensional arguments then indicate that the decay rate is of 
the form 

r ~ g?gm\ , (1.4) 

which is of order 10 -36 GeV. 

Despite this, we have two motivations for performing this calculation. First, in theories that 
extend the standard model of electroweak interactions and/or the gravitational interactions, 
including theories of Lorentz invariance violation, the rate might be different. This calculation 
can pave the way and serve as a test case for similar calculations in the context of those extended 
theories. Second, the methods employed here can also be helpful in the calculations of related 
processes where the Z boson appears as a virtual particle, such as vv — ► jQ, and it is conceivable 
that they can have physical relevance in some astrophysical contexts. 

The amplitude for the process is determined by a set of one-loop diagrams in perturbation 
theory, that we divide into two classes. The first class, to which we refer as the fermion loops, 
consists of the diagrams that contain fermion lines circulating in the loop. The second class con- 
sists of the diagrams in which the W bosons circulate in the loops, and in principle the diagrams 
that contain their corresponding unphysical Higgs bosons and Fadeev-Popov ghosts. Among the 
various vertices required to compute the diagrams, we need the gravitational couplings of the 
W boson as well as those of the Z and the photon. We adopt the point of view that, for each 
gauge boson V = 7, W, Z, those couplings are determined by the interaction Lagrangian 

= -Kh^TjW ? (L5) 

where is the expression for the energy-momentum tensor that is obtained from the linear 
expansion of the gravitational Lagrangian term 

J2f> = V=^°° , (1-6) 

with _5?o being the canonical expression for the bilinear part of the Lagrangian. Thus, for 
example, for the photon, 

^ i] = -\f^F^, (1.7) 

while for the W and Z it contains the corresponding mass term. This amounts to adopting 
the unitary gauge for the W and Z gauge bosons and therefore, for consistency, we employ 
the unitary gauge throughout. The amplitude that is obtained in this way has the following 
properties, which justify this approach: (1) it satisfies the transversality condition required by 
the conservation of the electromagnetic current, (2) it satisfies the analogous condition required 
by the conservation of the energy-momentum tensor, (3) the diagrams yield a finite contribution 
to the amplitude. In the presentation that follows, we consider the fermion loops first which are 
simpler, and therefore they allow us to introduce a set of techniques that are useful for treating 
the more complicated W loops. 
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The paper is organized as follows. In Sec. |^1 we show that Lorentz invariance together 
with the electromagnetic and gravitational transversality conditions imply that the on-shell 
amplitude can be expressed in terms of three form factors, two of which are non-zero only if the 
CP symmetry is broken at some level. The remainder of the paper is devoted to the calculation 
of the only form factor that is non-zero at the one-loop level. In Sec. EJ we enumerate the 
fermion loop diagrams and the corresponding Feynman rules, and we verify that they form 
a gauge invariant set in the sense that their total contribution to the amplitude satisfies the 
electromagnetic and gravitational transversality conditions. We then proceed to calculate their 
contribution to the form factors and, as expected, the two CP-odd form factors vanish and only 
the one that is CP even survives and it is finite. In Sec.0]we carry out a similar procedure for the 
W loop diagrams, with analogous results. As already mentioned, the calculations are carried out 
employing the unitary gauge for the W, Z propagators, together with the gravitational couplings 
that follow from the canonical expressions for the energy-momentum tensor. Finally, in Sec. 
we use the results obtained in the previous sections to compute the decay rate. In the course of 
the calculations we have used several Ward-like identities that relate the various gravitational 
vertices and other algebraic manipulations, we which have summarized in the appendices. 

2 General form of the amplitude 

We introduce the off-shell vertex function F\ pfiV (q,k), which is defined such that the on-shell 
amplitude for the process is given by 



tors for the photon and Z, respectively. Since we are calculating the on-shell amplitude, the 
momentum vectors satisfy the on-shell conditions 









0, 
Ml 





and they are related by momentum conservation 



p = k + q 




which in turn imply the kinematic relation 



2k-q = Mz- 



(2.6) 



The polarization vectors for the photon and Z satisfy 



e v [k)k v = 0. 





The polarization tensor for the graviton satisfies the analogous relations 
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and in addition it is symmetric and traceless, i.e., 



£ xp = S px , (2.10) 
£ Xp V\ P = 0. (2.11) 

Eq. (|2~TU|l implies that F Xpp v can be denned such that 

Fxppu = Fp\pi/ • (2-12) 

In addition, the conservation of the electromagnetic current and the energy-momentum tensor, 
which are consequences of the electromagnetic and gravitational gauge invariance, imply some 
additional properties of F\ ppv which will be useful in the explicit calculation of the amplitude. 
In the rest of this section, we explore the consequences of these conditions. 

The fact that the Z and the graviton are electrically neutral has two implications 1 : (i) the 
conservation of the electromagnetic current yields the condition 

k v F Xpill/ = 0, (2.13) 

(ii) the absence of tree-level diagrams implies that F Xppv can be expanded around k = 0. We 
exploit these properties by writing 

F\pp,v = ^XpfM, + & ^Xppua ' (2-14) 

where TP is independent of k. Since Eq. (|2.13j) must be satisfied for all k, it implies that 

^Txppu = 0, (2.15) 



' Xppu 

k v k a T x \ pua = 0, (2.16) 



which in turn imply 



Txpp.ua = antisymmetric in v <-> a . (2-17) 
Therefore, the amplitude can be written as 

M = E x P*ey va *t Xpilvol {q, k) , (2.18) 

where we have defined 

fua = k u E a k a E v . (2.19) 

and t\ pixva is some undetermined tensor. 

In analogous fashion, the conservation of the energy-momentum tensor yields the condition 

q x 3x P = o, 

q p jx P = 0, (2.20) 

1 When there are charged particles, Eqs. 12.131 and 12.141 are not valid. A relation analogous to Eq. 12.131 
holds when the charged particles are on-shell, but in that case the tree-level contributions (Born diagrams) render 
the amplitude singular at k — 0. 
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where we have defined 

jx P = e u *e»F XpfMJ (q,k). (2.21) 
Mimicking the argument of the electromagnetic case, we write 

hp = 3% + 3%X + jlparf? » ( 2 - 22 ) 

where j® and j\ are independent of q. As before, it follows that the transversality condition 
requires 

3x P = 0, 

3\pa = ~3ap\ ■ (2.23) 

In addition, Eq. (|2.12l) implies that 

JpXa = Jlpa • (2-24) 

By applying Eqs. (|2.23l) and (|2.24l) successively and repeatedly, we arrive at 

Jxpa = —jxpa > (2.25) 

or in other words j\ = 0. 

Considering now the structure of j\ paT , we can see that it has the following symmetry 
properties: 

1. It must be symmetric in the indices A, p because of Eq. (|2.12l) . 

2. It must be symmetric in the indices a, r because of its definition in Eq. 1)2.22(1 . 

3. Because of Eq. (|2.2U|) . it must be antisymmetric when we interchange anyone of the indices 
(A,p) with anyone of the indices (r, a). 

Moreover, any term in j\ par containing q\ or q p vanishes in the amplitude because of the on- 
shell transversality conditions of Eq. (|2.9() . and any term containing q a or q T vanishes as well 
in Eq. (|2.22j) since q 2 = 0. Thus none of the Lorentz indices A, /?, r, a can be attached to the 
momentum q. Remembering Eq. (|2.18jl . we know that j\ pTa must have a factor of f a p, and of 
course it also involves a factor of eft. Thus, combining these pieces of information, and excluding 
for the moment possible terms involving the Levi-Civita tensor, it follows that j\ pTa must be of 
the form 

fxpra = Fr XT Z pa + (\~p), (2.26) 
where F is some undetermined scalar, and we have introduced the notation 

Zpa = kp^cr ~ k<j£p > (2.27) 

or equivalently 

Z pa = Vp el- Va e z p . (2.28) 
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The terms involving the Levi-Civita tensor are most easily enumerated by noticing that all the 
requirements stated above for j\ pTa are satisfied if f a p or Z a p in Eq. (|2,26|) is replaced by its 
dual 



fa/3 
Z a f3 



c^aflpvf^ 



2^af3puZ IM , 



respectively. Thus the most general form for the amplitude is 



£ xp *a a t 



Ffl T Z pa + F x fl T Z pa + F 2 f{ T Z pa + (A <- p) , 



(2.29) 



(2.30) 



where F\ and F2 are two additional Lorentz scalars. Notice also that a term involving both / 
and Z is not included, since the product of two epsilon tensors can be written without it, and 
all such terms are already exhausted in Eq. (|2.30[) . 

It is convenient to write the expression for F\ piiv that follows from the form of the amplitude 
given in Eq. ([2.30)1 . Using the definition of Eq. i|2.1j) we get 



F\ PiW = (F(k\q v - k ■ qr] u \)(k p q pj - k ■ qr}p P ) 

+Fi[qk] Xv (k p q p - k ■ qr] pp ) + F 2 {k\q v - k ■ qrt v \)[qk] pi }j + (A <->■ p) , (2.31) 

where we have used the shorthand notation 

W = £^ a pq a k? . (2.32) 



3 Diagrams with fermion loops 

3.1 Diagrams and Feynman rules 

As already mentioned in the Introduction, we consider first the fermion loop diagrams, which are 
shown in Fig. ^ The Feynman rules for the various vertices that appear there have been given 
in the literature [71l51l9* HTUllll| . and for convenience they are summarized in Fig.[5J Any fermion 
/ that circulates in the loop must be electrically charged since it is attached to the photon line, 
and we denote its charge by eQf where e is the charge of the positron. The gravitational vertex 
function of the the fermion is given by 
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which can also be written as 
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V\ p (p,p') = - Jx(p + p')p + J P (p + p')x - -V\ P i> + i>' - 2m/ , (3.1 



where 



V\ P {p,p') = -\ax pa p(p + p') a l P + mfV\ P , (3-2) 



a\pnu = VXpViiu ~ 7}Vx P \pu ■ (3.3) 
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(e) (f) 

Figure 1: 1-loop diagrams for the process Z — > 7 + Q involving fermions in the loop. In the 
external lines, the thick saw-tooth lines refer to the Z-boson, the thin wavy lines to the photon, 
and the braided lines to the graviton. 

In writing this form, we have used the shorthand notation 

Vxp^u = VXfiVpu + rixvVp/i ■ (3-4) 

The tensor defined in Eq. Q3.3|) also appears in the coupling of a fermion bilinear to a gauge 
boson and the graviton. For example, if we denote Feynman rule for the Z-boson coupling to a 
fermion by — igj^/ (2 cos 6\y), where 

j^^Xf + Yfjs), (3.5) 

the fermion- Z-Q vertex is given by — ingaxpy^ / {2 cos Ow), as indicated in Fig. [21 In the 
standard model, 

Xf = T Lf - 2Q f sm 2 9 w , 

Y f = -T Lf , (3.6) 

where T^f is the eigenvalue of the diagonal generator of SU(2) acting on the left-chiral component 
of the fermion. 

As emphasized earlier, we employ the unitary gauge for the gauge bosons, which means that 
their propagators and the gravitational couplings are determined using the canonical form of 
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l \p 



-iKV Xp {p,p') 



h 



l Xp 



A a (k) 




Ap(k>) Z a (k) 




Z p {k') 



-iKC Xpa g(k,k') 



-^7/x/(2 cos W ) 



i\p 




-iKeQaxppul 11 



Z„ 




-ingaxppvY / {1 cos 9 W ) 



Figure 2: Notations for Feynman rules for couplings that appear in Fig. ^ The charge of the 
fermion is eQf. Various symbols appearing in this figure have been explained in the text. 

the kinetic term in the bilinear part of the Lagrangian. The gravitational vertex function c Xpa p 
of the photon is then given by 

c\p a/ 3(k, k') = r/xpiiiapk ■ k' - k' a kp) - r) a p(k\k' p + k' x k p ) 
+k(i{ri\ a k' p + r] pa k' x ) + k' a (r]xpk p + rjpph) 
-k • k'(r] Xa r] p p + r/ A(3 ?? pa ) , (3.7) 

and it is is useful to note that it satisfies 

k a c Xpal3 (k, k') = k'P CXpaf3 (k, k') = 0. (3.8) 

For our present purposes it is convenient to write it in the more compact form 

c Xpa p(k, k') = rj Xp (r] af3 k ■ k' - k' a k p ) - (j] Xp \ a p^ v - V\p\au,(3p) ^k w , (3.9) 

where we have introduced the notation 

flXp^P^u = r ]\p\af3 r ]p,v + r \\p\pv r la(i > (3.10) 

with r\ Xp \ pu defined in Eq. (|3.4|) , Similarly, the gravitational vertex function of the Z-boson is 
given by 



c \ pa p( k > k ') = C Xpal 3{k, k') - M z a' Xpa(3 



(3.11) 



where 



l \ppu ~ "Hxpilpu "Hxplp,!/ 



(3.12) 
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With this choice of vertices, the Z-propagator to be used is then given by 



Df(k) 



1 



k 2 - M 2 



k a k^ 



(3.13) 



For the photon propagator D a ^(k), the equation of motion that follows from the Lagrangian 
given in Eq. (|1.7|) determines only the transverse part of the solution, leaving an undetermined 
longitudinal part. That is, in the presence of a source J p , the vector potential is given by 



where 



k a k^ 



A: 2 



+ D L k a kP, 



(3.14) 
(3.15) 



with Dl being an undetermined scalar function of k. The consistency of the equation of motion 
requires that J a be conserved, which means that Dl drops out in Eq. (|3.14|) . and in fact that 
we can set 



D a ?(k) 



-T) 



a/3 



k 2 



(3.16) 



in the solution. The counterpart to this result in the context of our calculation is that, since the 
photon propagator enters only in the diagram in which the graviton is attached to the external 
photon line, then by virtue Eq. (|3.8|) for practical purposes the propagator can be taken as given 
in Eq. (|3.16|) . which is what we adopt. 

3.2 The amplitude 



( f ) 

We denote by the contribution to F\ ppv from one fermion / in the loop. Then writing it 

in the form 



7 (f) 

Xpfiu 



2 cos 



Ke 9 n T (f) 

^f ± Xpp,v 



(3.17) 



i if) 
the contributions from the various diagrams to are given by 

'%S(l - k) lv S{l)V Xp {l + q, l)S(l + q) 
%S(l)V Xp (l + q, l)S(l + q) lu S(l + p) 



j,(/:a) 
Xppv 


1 J 4 ( m 

= V(2.r Tr 


T (f-b) 
Xppv 


= 'Iw* 


T (f-c) 
Xppu 


= ax pai ,Il p a {p) 


T (f--d) 
Xppv 


= a Xppa^- a u(k) 


r (/ = e) 
Xppv 


= c Xpau (p,k)D 


T (f-i) 
Xppv 


= c z Xppa (p,k)D 



where 



n^(fc) 



(2vr) 4 



Tr 



%S(l)j v S(l + k) 



(3.18) 



(3.19) 
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Hu U (k), which is immediately recognized to be the fermion contribution to the polarization 
mixing tensor, satisfies the transversality condition 



Fn^(fc) = 0, 



(3.20) 



which is a consequence of electromagnetic gauge invariance and which can be proved explicitly 
by means of the elementary identity 



k v S{l) lu S(l + k) = S{1) - S{1 + k) 



This property of EL^fe) implies that it is of the form 

k p k u 



IV, (fc) = 77^ 



n(*), 



(3.21) 



(3.22) 



where H(k) is given by a (logarithmically divergent) integral that can be obtained from Eq. 
()3.19j) but whose precise value, as we will see, is not relevant for our calculation. 



3.3 The electromagnetic transversality condition 
Eq. (EHUj) implies 



k u T 



(/:d) 



Xppu 



ft 1 Xppv — u ) 



(3.23) 



and using Eq. ()3.8|) 



k v T {f:e) 



0. 



L xpixv " • (3-24) 
Then using Eq. Q3.21JI and remembering the momentum conservation equation Eq. ()2.5|) . we 
obtain 

d 4 ; 



ft 1 Xppv 



ft 1 Xppv 



(2vr) 4 

(2vr) 4 



Tr 
Tr 



%{S(l - k) - S(l)}v Xp (l, 1 + q)S(l + q) 
%S(l)V Xp (l, I + q){S(l + q)- S(l+p)} 



(3.25) 



and therefore 

ft 1 Xppv — 1 



d 4 l 

(2tt)< 



Tr 



%S(l - k)V Xp {l + q, l)S(l + q)- %S(l)V Xp (l + q, l)S(l + p) .(3.26) 



Changing the dummy loop momentum in the first term this can be written as 

d 4 Z 



h.u T (f : a+b) _ . 
K 1 Xppu — 1 



Tr 



(2^) 4 

and using the identity given in Eq. \KA 



%S(l){Vxp(l+p,l + k) - V Xp (l + q,l)}S(l+p) 



(3.27) 



uurp{f:a+b) 
" ± Xppv 



a X pauk u U p a (p) 



* ± Xppv ■ 



Using Eqs. (l3~2^|) . (l3~2i]) and ©29 > it follows that 



k u T U) 



Xppu 



0. 



(3.28) 



(3.29) 



which establishes the electromagnetic gauge invariance of this set of diagrams, in the sense that 
their contribution to the amplitude satisfies the requirement due to the conservation of the 
electromagnetic current. 
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3.4 The gravitational transversality condition 

In similar fashion, we now establish that 

gV*(fc)4(p)F^ = 0. (3.30) 

In contrast with the electromagnetic case, here both the Z and the photon are assumed to be 
on-shell. 

Using the identity 

S(l)q X V Xp (l + q,l)S(l + q) = l p S(l) - (l p + q p )S(l + q) 

+ \S(1) ( lp i - i lp ) - \ ( lp i - Hp) S(l + q) , (3.31) 

together with the Dirac-matrix identity 

lalplp + Iplpla = ^{jlaPlp + Tlppla ~ Vaplp) , (3-32) 

and making judicious shifts in the loop momenta in various terms we obtain 

q X Txp^ +h) = k p U^(k)-p p n^{p) + l(q v U w {p)-r, up q a U^(p^ 

+~ (q p n pu (k) - i] pp q a U au (k)) . (3.33) 
The definition of a\ ppv given in Eq. (|3.3I) allows us to write 



<fTxp£ = Qp^Ap) ~ ^quU pp (p) - -ri up q a n pa (p) , (3.34) 

q X Txp^ ] = q P ^u(k)-^Ii p u(k)-^ W a U au (k), (3.35) 

and therefore 

T Xp'^ +b+C+d) = PpU^-kpU^-v^U^-^Il^ik). (3.36) 
For diagram Fig. using Eq. (|2.5jl we first find that 

q X e u *(k)e^(p)c Xpa u(p, k) = e v *(k)e fl (p) (k a r) p „ - k p r} av )p 2 + (q u k p - r} pv k ■ q)p a , (3.37) 

where we have used Eq. Q2.8|) and the on-shell photon condition k 2 = 0. By the transversality 
property of IL^ we can use 

p a D^{p)U^(p) oc Pa U a p (p) = , (3.38) 

which implies that 

«V^A04(P)Z&? = e^(k)e^p)[k p U^p)- Vup k a U pa (p)] . (3.39) 
Similarly, for the expression involving diagram of Fig. we find 

q X e v *(k)e%(j>)c z Xppa {p, k) = e v *(k)e%(j>) [(p p r] pa - p a ri pp ){k 2 - M 2 Z ) + (q p p p - rj pp p ■ q)k a \ , 

(3.40) 
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and using the Z-propagator from Eq. (|3. 13|) 



which gives 



e^(k)e^p)cf pfWl ip,k)Df (k) = e v *(k)e»(p) 



riwQ 13 ~ VpP P 



(3.41) 



q X e»*{k)e%{p)T^ = e»*{k)e%(p) [^^{k) - Pp Il^(k)} . (3.42) 
Adding Eqs. (|3~3l?j) . (I3~3H1) and (033), and using Eqs. and (|3~277|) . Eq. (EOUl) is established. 



3.5 Calculation of the form factors 

The fact that the contribution of the fermion loop diagrams to the amplitude satisfies the elec- 
tromagnetic and gravitational transversality conditions implies that it must have the structure 
given in Eq. (|2.31|) . In order to extract the corresponding contribution to the form factors we 
need not evaluate in full the integral expressions for the various diagrams, and instead we can 
proceed as follows. The form factors F\ % are easily identified as the coefficients of the terms 
containing the Levi-Civita tensor. As we will see, no such terms appear so that these form 
factors are zero. For the form factor F, we can fix our attention on the contributions to just 
one of the terms that appear in Eq. (|2.31|) . To be specific, we choose Fk\k p q p q u . Since, by Eq. 

inzn>> 



(3.43) 



such terms will arise from the terms in the integrals that contain a factor of either k\k p q p q v or 
k\k p k^q u . A term of either form will be called a kkqq term for the sake of brevity, and these are 
the only ones that we need to track in the evaluation of the integrals. 

We now consider the contribution from each diagram to that kind of term, using the expres- 
sions given in Eq. (|3.18|) . Since T^^j and T^'J^ depend only on p and k, respectively, neither 
one contains a kkqq term. Using Eqs. (|3.13l) and ()3.22[) it is immediately realized that T^ 'ij and 



j-Xpuu are proportional to c\ ppu (p, k) and cf ppu (p, k), respectively, neither one of which contains 
a kkqq term as can be seen simply by looking at their definitions given in Eqs. (|3.9j) and (|3.11|) . 
Thus, we are left with t[^;^} and r}{,',? which using the graviton vertex in the form given 



XppU 



in Eq. (|3.2j) and remembering Eqs 



T. 



(/:a) 



Xppv 



P (/:b) 
" Xppv 



d*l 

(2vr) 4 
(2vr) 4 



IpTr 



Xppv 

and ()2.11j) . can be written as 

ft + m f)lvQ/ + m fhx(f +i + mf) 



[{I - kf 



m. 



+ ■ 



m 2 ](l 2 



m. 



IpTr 




j + m/)7A(f + mfhu(y+ ft + mf) 


[( 


l + kf 


- m 2 ][(l — q) 2 — m 2 ](l 2 — m 2 ) 



Changing the integration variable from I to —I in the integral for T^^y, 
cyclic property of the trace, together with the relations 



(3.44) 



and then using the 



C- 1 S T (£)C 



~1p 
S(-£) , 



(3.45) 
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and similar ones, it follows that the 7^75 axial current coupling from T^** and T\{.*„, are 



opposite, while the vector current term 7^ is the same. Therefore, 

d 4 / fx P ^u{l) 



' Xppu 



r (/:a+b) 
" \pfiu 



2iX 



where 



fx Pt w{l) = 1 p Ti ' 



(2iry [(I - k) 2 - m 2 ][(l + q) 2 - m 2 ](l 2 - m 



7 M (f- jk + m f )^(f+ m / ) 7A (f + i + m f ) 



(3.46) 



(3.47) 



The absence of the 75 term implies that there are no terms containing the Levi-Civita tensor, 
so that there are no contributions to the form factors iq,2- Continuing to D dimensions and 
parameterizing the integrals in the standard fashion we obtain 



T 



(/:a+b) 



Xppu 



AiXi 



dx 



l-x 



dy 



fa(l + xk-yq) 



m 2 + xyM\ 



(3.48) 



where we have used the on-shell relations given in Eqs. (|2.2|) . (|2.3|) and (|2.6|) . Evaluating the 
trace and focusing on the kkqq terms as described above, we get 



p (/:a+b) 
" \pp,v 



QUXfq^kxkf, 



(2tt) 4 Jo 



dx 







i-x x 2 y(l 
dy - 



x 



y) 



l 2 



m 



+ xyMl 



(3.49) 



where we have continued back to four dimensions since the resulting integral is convergent 
and the ellipses indicate that have omitted all the other terms. Performing the momentum 
integration and recalling the overall factors in Eq. (|3.17j) . we finally find the contribution from 
any fermion in the loop to the form factor F to be 



F U) 



Keg 



2tt 2 cos 9\v 



Q f X f I{m f /M z ) 



where 



1(A) 



d:X 







l-x x 2 y(l 
dy 



y) 



A 2 



xy 



(3.50) 
(3.51) 



4 Diagrams with W loops 
4.1 Diagrams and Feynman rules 

The one-loop diagrams involving the VF-boson are shown in Fig. |31 and the relevant Feynman 
rules are summarized in Fig. 0J As stated in the Introduction, we employ the unitary gauge, so 
that the propagator and gravitational vertex function of the W are given by 

c\pap(k> k ') = c\pap(k,k') - M^a' Xpa/3 , (4.2) 
in analogy with the corresponding quantities for the Z boson. The trilinear boson couplings 



13 




Figure 3: 1-loop diagrams for the process Z — ► 7 + Q involving charged gauge bosons in the 
loop. In the external lines, the thick saw-tooth lines refer to the Z-boson, the thin wavy line to 
the photon, and the braided lines to the graviton. 
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h\ p (q) 



IV+(p)vwwmWj(p) W+ipx) 




A u (k) A u (k) h Xp (q) 



-ieN a(3u {p 1 ,p 2 ) 



-iKeNx pa p v (pi,P2, ~k) 



W+ (pi ) wv/x&x/w (p 2 ) 



Zpkp) h \p{1) 




w;{p 2 ) 



-ig cos OwNapft (pi , p 2 ) 



-my cos 9 w N Xpa0ll {pi ,p 2 ,p) 




W+{P2) 



w + 




-ieg cos OwRap^u 



-IKCg COS 6 W R\pa(3iiv 



Figure 4: Notations for Feynman rules involving the VF-bosons in the unitary gauge and their 
extensions to include linearized graviton couplings. For the free VF-line, the notation represents 
the inverse propagator. 



involving the W have the following form in the momentum space, 

J^cubic = - 9W a W^N a ^( Pl ,p 2 ,p 3 ) , 



(4.3) 



where W a is the field which annihilates the W + boson, and W9 is the field operator for the 
neutral SU(2) gauge boson. The momenta have been written in the same order as the gauge 
bosons fields, with p 2 flowing out of the vertex and the other two flowing in, and 

N a Py(pi,P2,P3) = VP~tVaa{P2 + P3) a + V~farif3cr{pi - P3) a - VapV-yviPl + P2) a ■ (4.4) 

Since only two of the three momenta are independent in this vertex, we will often use the shorter 
notation 

N a /3j(pi,P2) = N a pj(pi,p2,P2-pi)- (4.5) 
Similarly, the quartic WW'jZ interaction in the flat space Lagrangian is given by 



^quartic = -eg cos9 w W ' a W^Z 'nAyRapnv 



where 



(4.6) 
(4.7) 
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Following the linear approximation to the gravitational interactions already outlined, we then 
obtain the corresponding vertices involving the graviton, characterized by the vertex functions 

N\ p af3j(Pl,P2,P3) = T]\ p N a p J (p 1 ,P2,P3) ~ VXp\f3-y,aa{P2 + PzY 

~ VXp\~/a,f3a(Pl ~ PzY + VXp\af3,~/a{Pl + PlY , (4.8) 

(4.9) 

as indicated in Fig. 21 where the symbol r/ with six indices was defined in Eq. ()3.1(J|) . 
4.2 The amplitude 

Using the Feynman rules just discussed, and defining T^^ v by 

^ = ^cos%rW, (4.10) 
the contributions to T^J U from the diagrams in Fig. 01 are given by 

«S£ a) = /(04^(m+pW(o 

x cY paT (l + q,l)D$(l + q )N 1&u {l+p,l + q)D${l + p) , 
ff££ b) = J^Na P p(l-k,l + q)D^(l-k) 

x N aru (l, I - k)D$(l)cY n5 (l + q, l)D$(l + q) , 
iTZu C) = J^N a ^(U+p)D%(l)N Xparu (l+pJ-k)D^(l+p), 

x J ^N a/3fl (l,l+p)D^(l)N aT1 (l+p,l)D^(l+p) , 
V$£P = 4 piia (p,k)Df(k) 

x J 4JjiN^p(l, I + k)D%(l)N aTV (l + k, l)D$(l + k) , 

iTiZl S) = Ra^J^D^{l)cf paT {l + qi l)D^{l + q), 

■ T (W:h) R /" d 4 / 

L1 Xpp,v — n-Xpafipv J ^vr) 4 W^li 

iT { Z? = c Xp UP,k)D^(p)R aT ^J 1 0^D^(l), 

«$£ J) = cZ ppa (p,k)Df(k)R aTpu J^D%(l). (4.11) 
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4.3 The electromagnetic transversality condition 

First, using Eq. (|3.8|) it is immediately seen that 

k»T^ = k»T^=0. (4.12) 
For the rest, we use the Ward identities that have been given in Appendix Thus, using Eq. 

(ESI) 



WT^f = -k»R XpaPlxv (4.13) 



and comparing it with the expression for ^ in Eq. (|4.11|) we conclude that 

Fr (^.d +h)=0 _ (414) 

Similarly, applying Eq. (|A.6|) . shifting the momentum in one of the resulting terms and then 
applying Eq. ()A.9|) . we find 

k » T ^) = 0. (4.15) 
Applying Eq. (|A.6j) to T^f> and T^ b) we obtain 

+ / T^fiWU +p)D%(l)cY paT (l + q,l)(D^(l + q) - L$(l+p) 

^N a ^(lJ+p)D$(l)(cW yS (l+p,l + k) - cf^il + q,lj)D$(!+p) 

(27r )4 (^W*. l+P)~ N a ^(l ~k,l + qj)D<%{l)cf p(TT {l + q, l)D%(l + q) . 

(4.16) 

From this, and using Eqs. (|A.7|) and ()A.9|) . it follows that 

^C a+b+c+s)=0 > ( 4 - 17 ) 

which complemented by Eqs. ()4.14j) and H4.15|) establishes the property 

^411 = 0. (4.18) 

4.4 The gravitational transversality condition 

In similar fashion, here we establish that T^J V satisfies 

gV*(fc)4(p)T^2 = 0. (4.19) 
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d 4 l 

+ ' 



In the formulas that we obtain below for the contraction of q with the various amplitudes T^^\ 
we omit writing the polarization vectors of the photon and the Z to simplify the notation, but 
it should be understood that the relations are valid in general only when the contractions with 
the polarization vectors as indicated in Eq. 1)4.19(1 are taken, and the on-shell conditions for the 
photon and the Z are imposed. 

We start with the identity involving the W gravitational vertex and the propagator in the 
unitary gauge, 



D^\l + q)q x c^{l + q,l)D^{l) = IpDff' (l)-(! P + , 



)D^\l + q) 



+ Vpa (f'D^{l)+q^Dt{l + 
which resembles Eq. ()3.31|) for fermions. Using Eq. I|4.20|l . we obtain 

d 4 l 



(4.20) 



ia x T {W:g) 

l( l l \ppv 



R 



afipv 



l p D$(l)-{l p + q p )D%{l + q) 



a/3. 



(2vr) 4 

+r }pa q a D^{l + q)+ ^B^{1) 



(4.21) 



The first two terms cancel, as can be shown by making a change of variables in the second one. 
Using Eq. (|4,7|) the remaining two terms can be written as 

d 4 / 



ia x T (W:g) 

l 1 1 \ppu 



(2^) 4 

Next, using the definition in Eq. ()4.9|) . it is straight forward to write 
d 4 / 



(4.22) 



ia x T (W:h) 

L( i 1 \ppu 



(2tt)< 



qpRaf3pvD^(l) - 2(r] up qp + r] pp q u ^r] a pD^ (I) - ^r]p V q a D^ p {l) 



+2q a {rj up D^(l) + ^^,(0) + 2g M £>JJ(Z) + 2q v D%(l) 
Further, using Eq. ()3.37|) . we obtain 

d 4 / 



(4.23) 



ia x T iW:i) 



kpRarpi/ Vvp^ RarpP (.Qi/kp T/i/pk • q) — Rcr T pf3 I 

A' _ « 



(27T) 



Similarly, using Eq. 1)3.41)1 . we can write 

VppQ Rcrrau P pR<jr pv 



ia x T {W:i) 

%< l 1 \ppu 



d 4 Z 



mil)- (4.24) 



(4.25) 



and summing up Eqs. (|1~2^|) . (|4~2l?)) . ((01]) and (JOSJ), we obtain 



d 4 ^ 



- (lukp - Vupk ■ q)^Rarp/3 



(27T) 



4 ctp 



p- 



d 4 / 



(27T)' 



D%(1). (4.26) 



Regarding , by means of Eq. ()B.7|) we can write 



ia x T (W - c) 

l 1 1 \ppu 



-N af3ll (l,l+p)D^(l)D^(l+p) 



(M 



k p N aTU (l + p,l) 



+l p q S Raru8 + P p N aru (l + k,l)- r] pa q S N Stu (1 + fc, /) 
-rjprq 6 N aS u{l +P,l + q)~ VpuQ 6 N arS {l +p,l) , 



(4.27) 



18 



and Eq. (|3.37f) allows us to write the result for Q^T^' as the sum of the two terms, 



• X T (W:el) 

ia X T (W:e2) 
l 1 1 \ppu 



kp1] u k 7] pv 



d 4 Z 



■N aTS {l+p,l)N aPlt (l,l+p)D^(l)D^(l+p) , 



\<lvkp - i] U pk ■ q) 



pi f d A l 



l N a/3p (l,l+p)N aT , / (l+pJ)D^(l)D^(l+p) 
(quk p - r, up k ■ q)^R aT „ s J tS^W^) » ( 4 - 2 8) 



p 2 J (2tt) 4 ' 
5 <■ d^_ 

(2tt)< 



where we have applied Eq. (|A.12j) in the last step. If we consider the sum q^T^ u ° +e1 ^ , the term 



from q^T^^J 1 ^ containing k p cancels an identical term from Eq. (|4.27j) . The remaining term 

from q^T^'^ combines with the term proportional to r\ pv from Eq. (|4.27|) . and applying Eq. 
(lA~12l) they yield 



(4.29) 



Substituting the expression for R a paS given in Eq. (|4.7jl . the term proportional p^ vanishes due 
to Eq. (|2,7|) and the remainder cancels one of the terms of Eq. (|4.26|) . The other term of Eq. 
(j4.26|) is canceled by the q^T^^J 2 ^ contribution, and in this way we obtain 



■\rp(W:c+e+g+h+\+i) 
%( 1 1 Xpp,u 



d 4 l 



4 N aP ^l,l+p)D^(l)D^{l+p) 
+P p N aTU {l + k,l)- Vpaq S N Sru {l + k, I) 
-r]p T q s N a Su{l+P,l + q) 



For the contraction of we use Eq. p. 41(1 to write it in the form 



Ipq RarvS 



(4.30) 



n 



N aps (l, I + k)N arv {l + k, l)Dff(l)D% (l + k), (4.31) 



(2tt) 4 

and applying the identity of Eq. (|B.5|) 

d A l 



ia x T {W:d) 

% 1 1 \ppu 



N aTl/ (l + k,l)D^(l)D^(l + k) 



{P P ~ l P )N a ^(l,l + k) 



+lpN a pp{l - q,l - q + k) - kpN af3p (l, I + p) - r] pa q S N S/ 3 p (l -q,l + k) 



-ri p pq S N a5p (l,l +p) - r] pp q 5 N a ps(l,l + k) 
Adding these results and using Eq. (|B.3|) . 



(4.32) 



. X T (W:d+i) 
l 1 1 \pp,u 



dn 



N aTU (l + k,l)D^(l)D^(l + k) 



IpQ Raf3p,6 



(27r) 4 " ' aTUK " ' , "i'i"w vi" w 
kpN a pn(l, l+p) - i]p a q S N S /3 p {l -q,l + k) - r) pj3 q 5 N a& p(l, I + p) 



.(4.33) 
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(a) 



(b) 



Figure 5: One-loop diagrams for the 7Z polarization mixing tensor involving internal VF-boson 
lines. 



Finally, for TjL ' and TjL ' we use Eq. (|4.20|) once again. By redefining the integration 
variable in certain terms, the results can be written as 



io A T (W:a) 
t( i 1 Xppv 



ia x T iW:b) 



d 4 l 

(2vr) 4 

+N aPtl {l -q,l + k)N aTU {l + k, l)D^(l + k)(- lpDff(l) + ri T5 q a Df p (l))UM) 



N a/ 3 (i (l,l+p)N aTV {l + p,l + q)D$(l+p)(l p Dff(!) + Ti a& q T Df p {l 

Sp 



d 4 l 



N a ^{l,l+ P )N aru {l + k, l)Dff(l) 



(l p + k p )D^(l + k) 



(4.35) 



(lp+P P )D^(l+p) + rfttfDfJl+p) + rf^DfAl + k) 



Adding Eqs. (|4~3U|) . (jH3HJ, (jOljl and (j4~35j) . and applying the identities of Eqs. (|B.3j) and (|B~4 

we verify that all the terms cancel, which proves Eq. Q4.19JI . 



4.5 Calculation of the form factors 



As we have already argued, the virtue of having proved that T^ p ^ u satisfies the electromagnetic 
and gravitational transversality conditions is that we are now assured that it has the structure 
given in Eq. (|2.31j) . Furthermore, there cannot be any contribution to the factors F\ and F2 
since the Feynman rules do not involve the Levi-Civita tensor. Therefore, to calculate the 
contribution to the form factor F, we can just look at the kkqq terms, as we did in the fermion 



case. Fortunately, as we now show, only T 



(W:a) 



Xppv 



and T 



(W:b) 



XpfMV 



produce such terms, and they are 

in fact equal so we need to evaluate only one of them. 

Since R a pp5 is independent of momentum, the amplitude of T^ pt f^ does not depend on k, 
and consequently does not produce a kkqq term. The same argument applies to the amplitude 



^\ P pv ' w hich is independent of k as well as q. Next, consider the amplitudes 7\ 
Their sum can be written in the form 



T (W:f) T (W:j) 
Xppv ' Xppu 



cl piMX {p,k)Df{k)Yl%\k) 



(4.36) 



where IT^P denotes the contribution of VF-loops to the "yZ polarization mixing tensor, repre- 
sented by the diagrams shown in Fig. [5J This tensor, as can be readily verified explicitly, is 
transverse and therefore it has the same form given in Eq. (|3.22|) . We can then apply the same 



argument given in Section 1531 to dismiss the T$f '„ and TA'^J amplitudes in the fermion case 



Xpfiv 



Xppv 
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to conclude that the combination T^'^ + T^Y' does not give any kkqq term, and a similar 

argument holds for T^' e ^ + as well. 

Considering T^' d \ the g-dependence can come only from the factor Nx pa f3p{lJ + k,p) in 
the integrand after we substitute p = k + q. Since 



N\ P afa(l, I + k,k + q) = {r]\p\ m ,pa ~ Vxp\pp,a<r)(f + terms independent of q , (4.37) 



we obtain 



(W:d) 



Xppv 



JlXp\pa,/3a - VXp\f3p ,(\<T 



d 4 l 



D^(l)N aTU (l + k, l)L$(l + k) + 



(4.38) 



where the dots denote ^-independent terms and therefore do not contain any kkqq terms. The 
remaining integral in this equation yields a three-index tensor that depends only on k. Since 
k 2 = 0, the only contribution from the integral that yields a term with four powers of momenta 
in Tj^ u d ' contains the factor k a k^k v . Eq. ()2.8j) then implies that it does not contribute to 

the amplitude, and thus we conclude that ^ does not yield any kkqq term. The same 

conclusion is reached for Tj^' c ' as well, for which it is only necessary to note that 



where we have explicitly indicated the (p, k) dependence in order to state the argument. 



(W:c). 



(4.39) 



Thus, only T 



(W:a,) 



Xppv 



and T- 



(W:b) 



Xppv 



can give the required type of terms. Furthermore, by changing 



the integration variable from I to —l + k in the expression for given in Eq. Q4.11J1 and 

using the symmetry properties of the couplings involved, it follows that 



r (W : a) 
" Xppv 



T 



(W:b) 



Xppv 



(4.40) 



Therefore we need to evaluate only one of these two, and we choose T^Y ^\ Wo lirst. write the 



amplitude it in the form 



iT 



(W:b) 



Xppv 



where we have defined 



and 



(2vr) 4 



A w (l - k)A w (l)A w (l + q)S. 



Xppv 



(4.41) 



(4.42) 



Sxppu = N al3fJ _(l - k,l + q)N aru (l,l - k)c Xn5 (l + q,l) 

ri s ~ 



_ v *r + (J-kr(i-ky 



s 



(o) 

Xppv 



+ ^Xppv + S 



M 2 W 

(4) 
Xppv 



-1] + 



M 2 W 



+ S 



(6) 

Xppv ' 



^ , {l + qf{l + qV 



M 2 W 



(4.43) 



In the last step, we divide the terms in Sx P pv into four terms, each one labeled by an index that 
denotes the number of inverse powers of M\y that it contains. 
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We consider first the term which has an overall factor of l/Myy, 

s { Zu = ^(!-k) a (!-krri'(! + qf(! + qy 

1V1 W 

xiV Q/3/ ,(/ — k,l + q)N aTU (l, I - k)cf nS (l + q,l). (4.44) 

To reduce this term and the other ones, we use the following identities 

a a N a/3l (a,b) = Pp 7 (b - a) - Pp 7 (b) , 
brN a py(a, b) = P ya (b - a) - P ya (a) , 

a a bPN aj3l (a,b) = a ■ (b - a){b - a) 7 - (6 - a) 2 a 7 , (4.45) 

where 

P a p(l) = -l 2 Vap + Up ■ (4.46) 

The identities in Eq. (|4.45j) . which are similar to the relations given in Eq. (|A.5|) . follow directly 
from Eqs. ()4.4|) and (|4.5[) . Applying them to the particular combination that appears in Eq. 
(|4.44|) . we obtain 

- k) T N aru (l, l-k) = l-kk u - k 2 l u , (4.47) 

which shows that does not contribute to the amplitude after using the on-shell conditions 

for the photon given in Eqs. ()2.8)1 and (|2.2j) . 

S\p U v contains three terms, with an overall factor of 1/Myy- One of them vanishes due to 
Eq. 1)4.47(1 . For the others, we use the definition of Eq. (|4.2j) and the identity of Eq. 1)3.8)1 to 
obtain another useful identity 

l' T c% aT {l, I') = -M^ra' Xpur , (4.48) 



which allows us to write 



5 £U = T72- a 'xpjsNa^{l ~k,l + q)N aTU (l, I - k) 



1 



M 



w 



r, aS {l - k) a (l - kf + V aT l a l S ] (I + qf(l + qy , (4.49) 

where in the last step we have used the definition of from Eq. (|4.2|) . along with the 

identity of Eq. 1)3. 8|) . Applying the identities in Eq. 1)4.45)1 . and omitting any term that does 
not contribute to the amplitude due to the on-shell conditions for the photon and the Z, the 
expression for S^} reduces to 



S aIuv - T7T a 'xfns( l + 9) 



M 2 z (l - k) p P 6 v {l) + l d {M 2 zVafl + P afl (l - k))P?(l - k) 



(4.50) 



Using the identity 

Pa»(l ~ k)P?{l -k) = -(l- kfP^l - k) , (4.51) 
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and noticing that we can make the replacements 



Vxslp ~ rjpsh , 



(4.52) 



because the neglected terms give a vanishing contribution to the amplitude when the on-shell 
graviton condition is imposed, we obtain 



S 



(4) 
Xppu 



M 2 W 



l\lp(l k)pl v 



2Mi 



(I - kf 



+ 



(4.53) 



where the dots denote terms which cannot produce any kkqq term. This can be reduced further 
by writing 



(l-k) 2 = A^(l-k) + M^, 



(4.54) 



and remembering that S^} is to be substituted in Eq. Q4.41j) . The Ay^(l — k) term then cancels 
with the A\y(l — k) in Eq. (|4.41|) . and the resulting integral, which depends only on q, does not 
produce a kkqq type term. Thus the relevant part of S-$ is just 



S 



(4) 
Xpfiu 



4Mf 



l\l p (l - k) ^ + 



(4.55) 



S\p , which from Eq. (|4.43|) is given by, 



? (2) 
Xpfiu 



N^l — k,l + q)N aru (l, I - k)c^ piS {l + q, I) 



M 2 W 



rf T rf\l + qf(l + g) 7 + rj aT r]^n S + rf s rf^{l - k) a (l - k) T ] , (4.56) 
can be treated in similar fashion. By using Eqs. (|3,8|) . ()4.45|) and (|4.52|) . it can be reduced to 



S 



(2) 
Xpfiu 



(P T p{p) ~ P T p{l ~ kj)N aTU (l, I - k)(r&l p + rf p l x 
+N al3l ,(l -k,l + q) - P?(l - k)) (rf x l p + 4h) 

+ ^t{ P » ~ P 2( l + ?)) { P '(k) ~ P*(l))cZ, s (l + 3,0 + • ■ ■ , (4-57) 



which can be simplified further by neglecting the terms that do not contribute to the on-shell 
amplitude and focusing on the kkqq type terms. Omitting the rest of those details, we finally 



obtain, for the relevant part of S 



(2) 

Xpfiu ' 



5 



(2) 
Xpfiu 



6l\l p (l - k)fj,l p 



Adding Eqs. Q4.55JI and (|4.58|) we then obtain 



? (2+4) 
Xpfiu 



4 2 



Mf 



l\l p (l - k)Jv H , 



(4.58) 



(4.59) 
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and it should be remembered that do es not contribute to the kkqq terms, so that it need 

not be considered. Substituting Eq. (|4.59|) into Eq. (|4.41|) and parameterizing the integral in 
the standard way, we obtain the result from this part, 



p (W:b) 
' Xpfiu 



8i\2 



(2+4) 



Mf 



(2tt)D 



n 



" .hj h "" Al ■ !J " ] , ■ .... (4.60) 



P-M*+xyM* 



where now 



(4.61) 



Choosing the kkqq terms as we have indicated earlier, we obtain the contribution from this part 
to the kkqq term, 



T (W:b) 

Xppu 



8i 2 



(2+4) 



Mf 
M 2 W . 



kxkpq^qp 



(2vr) 4 Jo 



dx 



o 



1 x ^ x 2 y(l — x — y) 



P-M^ + xyM\ 



(4.62) 



where we have continued the integral to four dimensions since it is finite. 

? (o) 

Xppv 



We now work the term in Eq. (|4.43[) . which is given by 



? (0) 

\ppLV 



^VV^a/^G -k,l + q)N aTU (l,l- k)c^{l + q,l) 



(4.63) 



From the definition of the cubic couplings in Eq. (|4,4j) . and using Eq. (jH.Hjl . the part that can 
possibly contain kkqq terms is 



S 



(o) 

Xppu 



rg(l + k + 2q) a -2r} m {k + 



1rg{l + q), 



2^k d - V d u (l + k) a + 2rf a l v c Xn5 (l + q,l) + 



(4.64) 



where we have also made use of Eqs. ()2.8[) and (|3.43j) . In addition, we have noticed that the 
term proportional to in s does not produce a term quartic in the momenta, since the 
cubic couplings are linear in momenta. From the definition in Eq. (|3.9j) . and excluding all terms 
which either vanish on contraction with the polarization factors or do not contribute to the kkqq 
type terms because they do not have enough factors of uncontracted momenta, we find that we 
can substitute 



V-yshlp + vx-ylpi! + q)s + VXS^lp + (A <-> p) ■ 



The rest of the calculation is just straight forward algebra, and the result obtained is 



c(0) 
Xppu 



+ <l)p l v + QpQu + 8(k x l p + l\k p ) l^q u - q^k 



+ 



(4.65) 



(4.66) 



Substituting this expression into Eq. (|4.41|) we then obtain 
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(4.67) 



24 



which has to be added to the contribution obtained in Eq. 1)4.62(1 to obtain the complete ex- 
pression for the kkqq terms in T^^° '. Remembering that T^^f * and T^^^ give identical 
contributions, and recalling the overall factors defined in Eq. (|4.1U|) . we obtain the W contribu- 
tion to the form factor F, 

Fiw) _ ^^_^ 1{Mw/Mzh (4 , 8) 

where 1(A) has been defined in Eq. ((3.51(1 . 



5 Decay rate and discussion 

The on-shell amplitude is parameterized by the three form factors defined in Eq. 1)2.31(1 . Our 
calculations in Sec. Eland Sec. 0] show, first of all, that 

F l = F 2 = 0, (5.1) 

which is a consequence of CP invariance. The terms containing F\ and F2 in Eq. ()2.31() contain 
the Levi-Civita tensor, and are therefore odd under time reversal and under CP. In the standard 
electroweak model, CP violation enters any amplitude only through the charged current inter- 
actions of fermions, which do not appear in the one-loop amplitudes for the present process. 
Hence, at the one-loop level, the amplitude is CP conserving, and Eq. <|5.1|) reflects that fact. 

The only non- vanishing form factor in one-loop is F, for which the results given in Eqs. 
()3.5U() and 1)4.68(1 are combined to give 



neg 



4vr 2 M 2 cos 9w 



cos 2 e w U 1—] I(M W /M Z ) - 2j2QfXfI(m f /M z ) 

\ COS (/ \Y J j 



(5.2) 



where / is the integral defined in Eq. ()3.51() . The integral cannot be performed analytically, but 
we can make some approximations that are sufficient for our purposes. In the denominator of 
the integrand the combination xy has the maximum value | within the range of integration. We 
thus calculate the integral in two extreme cases, 

1(A) = 24 * (5 . 3) 

-7T for A > i. 

I 360A 2 4 

For all the fermions except the top quark, we use the first form, whereas for the top quark and 
the W in the loop, we use the second. Thus, 



F Ke 9 



cos 9w 



If 1 \ 5 / 2Ml\ 10 ~ / Ml 

6 ^— + — 1 % sin 2 6 W 1 ^ 

360 V cos 2 #W 12 I 15m 2 / 9 I 75m 2 



(5-4) 



where we have used the mass relation My/ = Mz cos 6w- Using e = gs'mOw and sin 2 9y/ = 0.23, 
this gives 

,2 



"z 



0.4-^—2 . (5.5) 
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The decay rate is determined straightforwardly from Eq. (|2,31|) , Using the familiar po- 
larization sum formulas for the photon and the as well as the corresponding one for the 
graviton|12j. 



(5.6) 



pol 




(5.7) 



pol 



(5.8) 



we find 



r = 



M 7 Z F 2 



(5.9) 



96vr 



and from Eq. (15.5(1 



r = 



o.i 



o 




(5.10) 



This result confirms our expectation in Eq. (|1.4|) about the smallness of the rate. However, 
as already mentioned in the Introduction, one of the primary motivations for performing the 
calculation was to understand some of the intricacies involved and resolve some of the technical 
complications in a way that can be used in similar, perhaps more complicated, calculations. 
In this sense, the proof of the general form of the amplitude and its parameterization given 
in Eq. (|2.31|) . together with the all the identities, manipulations and tricks that we have used 
both in the explicit calculation of the form factors as well as in the proofs of the consistency 
conditions (the electromagnetic and gravitational transversality conditions) are useful in their 
own right, independently of the fact that we have applied them in the particular context of the 
Z decay. It is particularly enlightening the fact that the use of the unitary gauge did not lead 
to any of the inconsistencies that are sometimes attributed to using that gauge. In fact, as 
we showed, the diagrams calculated with this gauge give an amplitude that is consistent with 
the transversality conditions implied by the electromagnetic and gravitational gauge invariance, 
which in turn allowed us to determine the amplitude and calculate the relevant form factor in a 
systematic and consistent fashion. By having considered this simpler system, it has allowed us 
to understand and develop some techniques that we believe can be useful for considering more 
complicated processes. 
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Appendices 



A Ward-Takahashi identities 

The simplest Ward-Takahashi (WT) identity familiar to us through QED, which relates the 
fermion two-point function to the fermion-photon vertex. Diagrammatically, it can be written 
as 



eQ 



I + r 



I 



( t rV \ 



) 




l + r I 



(A.l) 



where the fermion has charge eQ. The tree- level version of this equality was presented in Eq. 
1)3.21(1 . After the introduction of gravity, vertices involving gravitons appear in the theory. In a 
similar fashion, one can now prove the WT identity 



eQ 




\ 



l + t l+r 



\ 



( 



J 




l + s I 



\ 
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J 



i + r* i 



V 



(A.2) 



In this and other relations in this appendix, we use the shorthand 

t = r + s , 



(A.3) 



where r and s are the momenta of the photon and the graviton lines respectively, both considered 
outgoing. Recalling the definition of the vertices given in Fig. |5J we see that at tree level, the 
diagrammatic identity of Eq. (JA.2JI implies the relation 



V Xp (l + r + s,l + r)- V Xp (l + s, I) 



(A.4) 



which can be easily checked from the expressions for these vertices appearing in Eqs. (|3.1f) and 
(HE)). 

The general nature of these WT identities are now clear. For the charged W + bosons, we 
can write similar identities. Of course, these results depend on the gauge condition, and involves 
the unphysical scalar fields as well as ghost fields. In the unitary gauge, however, the unphysical 
scalar fields and the ghost fields are not present, and the relation look particularly simple. Here, 
we summarize some relations of this sort. In the diagrams which appear within the equations 
below, the horizontal lines denote W + , where the left line has incoming momentum and the 
right line has outgoing. The Z-line appears in a saw-tooth pattern, with an inward momentum 
p. The momentum convention for the photon and the graviton lines have already been stated. 

First, we can have the identity from a diagram similar to that in Eq. ()A.1|) . with the fermion 
lines replaced by the H^ + -lines. At the tree level, this will imply 



D Zp( l + r )~ D *W) = MafhO + r,l). 



af3 



(A.5) 



This can also be written as 



rW aa (l + r)N aPl {l + r, l)D^{l) = D aT {V) - D aT (I + r) . 



(A.6) 
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Similarly, replacing the fermion lines by T^ + -lines in Eq. (|A.2|> . we obtain another WT identity, 
which in the tree-level reads 



+ r + s,l + r)- cf M + s, I) 



r 7 iV 



Apa/37 



(I + r + s,l, -r) 



(A.7) 



We can have an extra Z-boson present in all diagrams involved in the identity of Eq. (|A.5|) . 
which will give us the relation 



'4~) - U2 

I ~ + t\ \l-rl + s 




(Ai 



At the tree level, this implies 

N a j3jQ, l + r + s)- N a p 7 (l -r,l + s) = - r 5 R a ^ s . 
We can also write down a similar WT identity involving the graviton, viz. 



(A.9) 




(A.10) 



The tree-level version of this identity reads 

Nxpap^U l + r,r + s)- N Xpa /3^(l -r,l,r + s)= - r l 'R\ pa /3 P v 



(A.ll) 



It is also instructive to see how some of these contractions behave under the integration over 
the loop momentum. For example, using Eq. HA.6|) . we can write 

(2vr) 4 



N arS (l, I + r)D^(l + r)N a ^(l + r, l)D&(l) 
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(2vr) 4 
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T RarSoi 



N utS {1, l + r)- N arS (l - r, l)\ D%(1) 

d 4 i 



(2vr)- 



(A.12) 



where we have used Eq. (|A.9|) in arriving at the last step. Through an exactly similar kind of 
argument, we can prove the relation 



d 4 l 

(2vr) 4 



NxparsV, l + r,r + s)D<%{1 + r)N a ^(l + r, I, -r)D$(l) 

= ~ r a R XpaT5a J 7^4^(0 • (A.13) 
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B Some other relations between various couplings 

In Appendix we considered the contraction of various vertices with the photon momentum. 
In this Appendix, we are considering contractions of a more general kind, in particular relation 
involving the graviton momentum. 

From the definitions of these vertices in Eqs. (|4.4I) and (|4.7|) . it is easy to see that one can 
write a relation between the cubic and the quartic gauge couplings: 

N a/ 3 7 (a, b) = a 5 R ia/ 3s + b 5 Rp ia s . (B.l) 

It is easy to see that Eq. (|A.9(I can be derived from this relation by making use of the identity 

Ra^S + Rf3-ya6 + R-ya(3S = , (B.2) 

which follows trivially from the expression for the quartic gauge coupling in Eq. (j4,7j) . The 
following identities, used at different stages of the calculation, can also be derived from Eq. 
(EH): 

iV Q/ 3 7 (a, b) - N al 3 7 (a, b - r) = r s Rp 7a s , (B.3) 
N a /3y(a,b) - N aj3l {a - r,b) = r 5 R loll35 . (B.4) 

Another important relation involves the contraction of the coupling N\ pa ^ with the graviton 
momentum. From the definition in Eq. (|4.8[) . it easily follows that 

q X Nx P af3~f {a,a + b,b + q) = (q p -a p + b p ^ N a ^ (a, a + b) 

+a p N a p 1 (a - q,a + b- q) - b p N a p 1 (a, a + b + q) 
-q S {v pa Nsi3y(a -q,a + b) + r] pl3 N a Sy(a, a + b + q) 

+Vp~fN al3 s{a, a + bfj . (B.5) 

Sometimes alternative forms of this identity is more useful, such as 

q X N Xpa p 1 {a, a + b,b + q) = q p N a ^ y (a, a + b) + q S (a p R a ^ s - b p Rp ia& ^j 

-q 5 {vp a Nsf3y(a -q,a + b) + i] p pN a&1 {a, a + b + q) 

+rj fn N a/3 s(a, a + b)), (B.6) 

or 

q X Nx pa pj(a, a + b,b + q) = (b p + q p )N af 3 y (a, a + b) - b p N a /3j(a, a + b + q) + q S a p R a p 1 8 

-q 5 (j]p a N Sl 3~,(a -q,a + b) + r] p pN a s^(a, a + b + q) 

+rt fr/ N a ps(a, a + b)), (B.7) 
which can be derived from Eq. ()B.5|) by using the identities of Eqs. ()B.3|) and (|B.4|) . 
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